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1. Introduction and Discussion 

The flow of non-relativistic fluids is described by the Navier-Stokes equations 

it + v. Vv = -VP + vV 2 v + f 

(1-1) 

V.v = 

where v is the fluid velocity, P the fluid pressure, v the shear viscosity and / an externally specified 
forcing function. Although these equations describe a wide variety of natural phenomena (see e.g. [1] 
) and have been intensively studied for almost two centuries, their extremely rich phenomenology 
remains very poorly understood. In particular, most fluid flows go turbulent at high Reynolds number, 
i.e. in the regime in which the viscous fluid term is negligible compared to the nonlinear convective term 
in (1.1). Although turbulent flows appear complicated and statistical in nature, it has been suggested 
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(see e.g. [2] ) that these flows are in fact governed by a new and simple universal mathematical structure 
analogous to a fixed point of the rcnormalization group flow equations. A completely new angle on 
fluid dynamics could well be needed in order to uncover such a structure. 

Recent investigations [3-13], within the framework of the AdS/CFT correspondence of string 
theory [14] have revealed an initially surprising relationship between the vacuum equations of Einstein 
gravity in an asymptotically locally AdS^+i space and the equations of hydrodynamics in d dimensions. 
More concretely these papers study a class of regular, long wavelength locally asymptotically AdS^+i 
solutions to the vacuum Einstein equations with a negative cosmological constant. These solutions 
are shown to be in one to one correspondence with solutions of the d dimensional hydrodynamical 
equations V M T M,y = 0. In the last equation the stress tensor T^" is a holographically determined 
functional of a d dimensional fluid velocity w M and temperature T. In the long distance limit under 
consideration it is appropriate to expand the stress tensor in a power series in the boundary derivatives 
of the velocity and temperature fields. Schematically 



where T is the local fluid temperature and T£ v is a local function of the fluid velocity and temperature 
of n th order in spacetime derivatives. The expressions for T£" for n < 2 have been explicitly deter- 
mined in the references cited above (see [13] for the most general result) and constitute a relativistic 
generalization of the incompressible Navier-Stokes stress tensor. In summary, classical asymptotically 
AdSd + i gravity is 'dual' to relativistic generalizations of the Navier-Stokes equations at long distance 
and time scales. 

Many theoretical and experimental investigations of fluid dynamics study the actual incompressible 
Navier-Stokes equations (1.1). It is consequently of interest to find a dual description of the Navier- 
Stokes equations (1.1) themselves rather than their their relativistic generalizations. This may be 
simply achieved by taking the appropriate limit of the results of [4-13], and this limit is the topic of 
our note. 

In order to find the gravitational description dual to (1.1), we adopt a two step procedure. First, 
purely at the level of fluid dynamics we make a straightforward and possibly well known observation. 
We note that the non-relativistic incompressible Navier-Stokes equations (1.1) are the precise and 
universal outcome of a particular combined scaling limit (one in which we scale to long distances, long 
times, low speeds and low amplitudes in a coordinated fashion) applied to any reasonable relativistic 
equations of hydrodynamics, i.e. the hydrodynamics a relativistic fluid with any reasonable equation 
of state 1 .The equations of fluid dynamics become non-relativistic at low speeds for the usual reason, 
and also become effectively incompressible, as we go to velocities much lower than the speed of sound 
(see for instance [1]). 

To be more specific we show that the equations of hydrodynamics reduce, in a precise fashion, to 



1 For instance the fluid that is dual to gravity studied below whose equation of state is dictated by conformal invariance. 
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(1.2) 
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the incompressible non-relativistic Navier-Stokes equations (1.1) under the limit 




(1.3) 



c 







where 5x is the spatial length scale, St the temporal scale while v % and SP represent estimates of 
the magnitude of velocity and pressure fluctuations about an ambient configuration of equilibrium 
fluid at rest. The rough contours of this choice of scaling are quite intuitive. It is clear we have to 
scale to long distances to be in the fluid dynamical regime. The dispersion relation for shear waves, 
cj = ivk 2 suggests that time intervals should scale like spatial intervals squared 2 . Scalings of distances 
and time intervals determine the scaling law for velocities. Finally the pressure variations are scaled 
appropriately to ensure that they cannot accelerate the fluid into velocities outside this scaling limit. 

As we have explained, the scaling (1.3) of the equations of relativistic fluid mechanics leads to 
the Navier-Stokes equations (1.1). It follows that this scaling operation is a symmetry of the same 
equations. This is easily directly verified. In particular, if the fields v l {x,t) and p(x,t) obey the 
unforced Navier-Stokes equations, then the rescaled fields ev l (ex,e 2 t) and e 2 p(ex,e 2 t) also obey the 
same equations. Consequently, the scaling operation described above is a symmetry of the unforced 
Navier-Stokes equations. 

As we have described above, the Navier-Stokes equations may be obtained as the scaling limit 
of any reasonable relativistic equations of fluid dynamics. In describing the connection with gravity, 
in most of the rest of this note, we will take the parent fluid dynamical theory to be conformal. It 
is natural to wonder how much of the full relativistic conformal group descends to a symmetry of 
the Navier-Stokes equations, and in what form it does so. It is obvious that the relativistic Poincare 
group descends to the Galilean symmetry group of the Navier-Stokes equations. In the next section 
we explain that the scaling symmetry operation described in the previous paragraph is loosely related 
to the dilatation operator of the parent relativistic conformal theory. 3 Further we demonstrate that 
all spatial special conformal transformations also descend to exact symmetries of the Navier-Stokes 
equations. After our scaling these transformations effectively turn out to be the 'boost' to a uniformly 
accelerated frame; the inertial forces one has to deal with when working in a non-inertial frame are 
compensated for by a shift of the pressure. These spatial special conformal transformations, together 
with the Galilean group and the scaling symmetry described above, form the (d + 2)(d + l)/2 — 1 
dimensional symmetry algebra of the Navier-Stokes equations. We list the commutation relations of 
this algebra 4 , and the action of its generators on the velocity fields, in detail in the next section. 

The conformal symmetry algebra described above is just subset of the full infinte dimcnsonal 
symmetry algebra of the Navier Stokes equations [16] (see also [17-19] for other related work) . The 
additional generators of the full symmetry algebra are very easy to describe; they consist of boosts to 

2 The scaling to arbitrarily low velocities projects out sound waves with dispersion relation u) <x k. 

3 This loose 'descent' is analogous to the relation of the dilatation operator of the Schrocdingcr group to the generator 
of scale transformations of the masslcss Klein Gordon equation. 

4 The nonrclativistic conformal symmetry group of the Navier-Stokes equations includes the contraction of spatial 
conformal generators Ki but does not include a generator that descends from the temporal conformal generator Ky. 
Our algebra is distinct from the Schrocdingcr group studied for example in [15] 
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a reference frame whose velocity is homogeneous in space but an arbitrary function of time. Just as in 
our discussion in the paragraph above, the pseudo force from such a frame change may be cancelled 
by an appropriate shift in pressure, and so is a symmetry of the Navier Stokes equations. 5 

The scaling limit described above admits an interesting generalization. Consider the equations of 
fluid dynamics on a base manifold G^ v = + where H^ v is small. By taking all terms that 
depend on iJ M „ to the RHS, V M T' il/ = reduces effectively to the equations fluid dynamics on the 
base space g M „ forced by an H^ v dependent forcing function. If we combine the scaling described in 
the paragraph above with the H 00 , Hij = 0(e 2 ) and H i0 = eAi(x l ,t), the effective resultant forcing 
function survives and is finite in the e — > limit. It turns out that the effective scaled forcing function 
depends only on Ai(x l ,t) and has a very simple form. It is precisely the force applied on a charged 
fluid by effectively a background electromagnetic potential A = 0, Ai = Ai(x l ,t). Consequently the 
'magnetohydrodynamicaP Navier-Stokes equations (i.e. the Navier-Stokes equations with a forcing 
function from an arbitrary background electromagnetic field) follows as a universal result of a scaling 
of the equations of relativistic hydrodynamics with small metric fluctuations. 

We now return to the duality between gravity and fluid dynamics. We apply the scaling limit 
described in the previous paragraphs to the equations of fluid dynamics that are holographically dual 
to gravity. This procedure gives us an asymptotically locally AdS^+i gravity dual to any solutions of 
the magnetohydrodynamical Navier-Stokes equations. The resultant metric describes small - but non 
linearly propagating fluctuations about a uniform black brane. We present the explicit form of the 
resultant bulk metric in section 4 below. 

It follows from our discussion that every solution to the Navier-Stokes equations (1.1) is also 
a scaling limit of a solution to Einstein's equations with a negative cosmological constant with one 
important caveat. The proviso is that many actual solutions of the Navier-Stokes equations describe 
fluids subject to hard wall type boundary conditions, and we do not (yet?) understand how to 
generate gravitational duals of these boundary conditions. However the qualitative effects of boundary 
conditions may easily be mimicked by appropriate forcing functions which are completely in our 
hands. In particular, while several experiments that study fully developed steady state turbulence do 
so in fluids with hard wall boundary conditions, there should be no barrier to setting up the same 
phenomenon for a fluid with no boundaries (e.g. on R d or on a compact manifold) with the appropriate 
forcing function. In order to make this expectation concrete we have identified one forcing function 
(of a likely infinite plethora of possibilities) applied to a fluid on R d , whose steady state end flow we 
expect to be turbulent at asymptotically high Reynolds numbers. In the rest of this introduction we 
describe this forcing function, and the reason we expect the flow it generates to be turbulent. 

Consider the forcing function A x = ae luJot+lkoV + cc, together with Ai — for i ^= x acting on (for 
concreteness) a fluid, on R 2 ' 1 whose spatial sections are parameterized by the Cartesian coordinates 
x and y. This gauge fields sets up a time and y dependent electric field in the x direction, together 
with a magnetic field in the plane. It is very easy to find one exact solution to the equations of fluid 
dynamics subject to this forcing function. On this solution v x is proportional to A x and v y = 0. 
This solutions describes a fluid driven into motion in the direction of an applied electric field (the x 
direction), while Lorentz forces from the magnetic field, in the y direction, are balanced by pressure 
gradients. All nonlinear terms in the Navier-Stokes equations vanish when evaluated on our solution. 
However nonlinear terms have an important effect on the dynamics of small fluctuations about our 
solution at high Reynolds number. In particular, in Appendix A below we demonstrate that nonlinear 

5 We thank J. Maldacena and R. Gopakumar for suggesting a symmetry enhancement along these lines, and thank 
J. Maldacena for drwaing our attention to [16], see especially the top of pg 68. R. Gopakumar and collaborators are 
currently further studying this algebra and its extensions. 
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terms drive some fluctuation modes (with momentum in the x direction) unstable at high enough 
Reynolds numbers. In order to explain the possible significance of this instability, it is useful to recall 
the usual situation with fluids at high Reynolds numbers. 

As we have described earlier in our introduction, fluid flows at high Reynolds numbers are very 
rich, and of potential theoretical interest. Nonetheless, in every highly symmetrical situation, there 
exists a simple 'laminar' solution, that preserves all the symmetries of the problem, at every value 
of the Reynolds number. This solution is often simple to determine analytically, and certainly shows 
none of the fascinating phenomenology of turbulence. However in interesting situations this solution 
becomes 'tachyonic' i.e. goes unstable to linear fluctuations that break some of the symmetries of the 
problem, above a critical Reynolds numbers. The 'end point' of this 'tachyon condensation' typically 
has richer dynamical behavior than the original solution itself. As the Reynolds number is further 
increased, further instabilities are usually triggered, and at arbitrarily high Reynolds numbers flow is 
turbulent. 

We suspect that the dynamical pattern described in the previous paragraph applies to the exact 
solution described in this paper. As we currently have no theoretical tools to predict the onset of 
turbulence in any fluid flow this is necessarily a guess, but one that we believe is natural, given the 
results of our stability analysis. This guess suggests that the stable steady state solution to AdS4 
gravity, with the effective gauge field A x described above, is dual to a turbulent fluid flow at high 
Reynolds numbers. Of course the particular situation described above is only one of a plethora of 
possibilities. We describe this solution in detail in Section 5 and Appendix A below only in order to 
have one concrete example of a gravitational set up that is likely to be dual to a turbulent fluid flow; 
not because we think that our particular is distinguished in any way. 

We believe it is likely that the fluid - gravity map will lead to interesting new insights on the nature 
of solutions of Einstein gravity in the presence of a horizon 6 . It also does not seem impossible that 
gravitational techniques and methods will prove useful in bringing new insights into the investigation 
of fascinating fluid phenomena like turbulence. In particular, it would be fruitful to understand the 
Kolmogorov laws on well-developed turbulence and their modification within the gravity framework. 
The symmetry algebra (3.5) may also throw new light on these issues. We leave further investigation 
of these issues to future work. 

Note Added : While we were completing this paper we received the preprint [20] which has 
substantial overlap with Section 2 of this paper. 



2. Scaling the Navier-Stokes Equation 

2.1 The Navier-Stoke Equation as a universal limit of fluid dynamics 

In this section we will display a scaling limit 7 that reduces fluid dynamical equations of the form 
V M T M,y = 0, to the incompressible non-relativistic Navier-Stokes equations usually studied in fluid 
dynamics text books (see e.g. [1]). 



"For instance, consider a black brane in asymptotically flat space. At sufficiently low energies, a beam of gravitons 
shot at this brane perturbs the non normalizablc boundary conditions of the effectively AdS near horizon region of the 
brane. It may be possible to choose this perturbation to drive the gravity in the near horizon AdS region turbulent. 

7 This section was worked out in collaboration with J. Maldacena. 
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The equations of relativistic fluid dynamics are 



v^r"" = o 

= pu^u v + pv^ - - cev^ + ... 

= g» v + uV, (2- 1 ) 
^ = V n* v »[} ^ VqM/3 + V/3Uq _ j^d.v^j , 6 = V% 

Here P is the pressure, p the energy density, r\ the shear viscosity, £ the bulk viscosity of the fluid, 
and g a p the metric of the space on which the fluid propagates. Each of the fluid quantities listed 
above may be regarded as a function of the fluid temperature - and also of chemical potentials if the 
fluid is charged. The ... in the equation above refer to terms of second or higher order in spacetime 
derivatives. 

If the fluid is charged, then we must supplement the equation (2.1) with an equation of charge 
conservation for every conserved charge. Below, we comment briefly on the scaling limit of this 
equation. 

2.2 A scaling limit 

Now let us study the motion of a fluid on a metric of the form G M „ = g^ + H^. We assume that the 
background metric has the form 

g^dx^dx" = -dt 2 + g lJ dx l dx : >. (2.2) 

(this is simply a choice of coordinate system, for a large class of metrics) while the small fluctuation 
H^ v is completely arbitrary. As we have explained in the introduction, we intend to view the fluid 
flow on the space with metric G M „ as an effectively forced flow on the space with metric g^ v . In order 
to do this we need to map the velocity field on the space to a velocity field on g^ u . This map, 
which must respect the requirement that u 2 = u 2 = — 1, may be chosen in a natural fashion if one 
thinks of the velocity field as generated by the path of 'particles' through spacetime. We now describe 
this in more detail. 

Let the fluid velocity on the space be given by where 



V 2 

where V is a d — 1 spatial vector with components V 1 , V a is a d component object with components 
(1, V) and V 2 is G a pV a V 13 (the indices a,(i run over d spacetime indices while the indices i,j run 
over the d— 1 spatial indices). Expanding to first order in H a p we have 

^ =1^ + 5^ + ... 

U " = y/l- 9ij VW ( 1,V ) (2-4) 
2 

where is the d velocity of the fluid referred to the metric g^. All terms in Su^ above depend on 
the fluctuation metric H a p; below we will take all these terms to the RHS of (1.1) and view them as 
contributions to the effective forcing function. 
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One obvious solution to the equations of fluid dynamics on the space with metric g^ v is simply 
the fluid at rest with constant pressure Po and density po- We now turn to a particular kind of small 
amplitude and long distance fluctuation about this uniform fluid at rest and with pressure Po and 
energy density po on a manifold that is 'close' to g^ v . More specifically we set 

Poo = € 2 h o(ex\e 2 t) 
H 0i = eAi(ex\e 2 t) 
Hij = e 2 h lJ {ex\e 2 t) 
V i = ev l (ex\e 2 t) 

P ~ Po 2 I i 2 , \ 

——=ep(ex ,e t) 
Pa + r'o 

and take e to be arbitrarily small. Although we have explicitly listed only the scaling of the pressure 
P above, the energy density p and viscosity v also scale in a similar fashion (this is consistent with 
the fact that they are all functions of the same underlying variables). We have normalized pressure 
fluctuations by p a + P rather than P for future convenience. 

We will now examine what happens to the Navier-Stokes equations under this scaling. Let us first 
start with the or temporal component of these equations. It is easy to check that 

V,T^ = e 2 [p e (vy)]+0(e*) 

Pe = PO + Po 

Consequently, in the limit of small e this equation reduces simply to V^w* = 0, where V, is the covariant 
derivative with respect to the purely spatial metric gij . 

Let us now turn to the spatial Navier-Stokes equations. After some calculation we find 



V M T^ = e 3 



Pe V l p + Pe V M (vW) - 2r ? V j I — J - CV, V.v - f 

2 \fg 



(2.7) 



where = (1, v 1 ). Using the equation V^u* = °^^ V - = 0, the coefficient of e 3 in the equation above 
may be simplified to 

V> + d Q v l + v.Vv 1 - v (V V + Ry) = ^£ _ d A l + F]v 3 (2.8) 



where = diAj — djAi is the field strength for the vector field Ai and v = r]/p e is the 'kinematical 
viscosity' of the fluid. 8 (2.8) takes a somewhat simpler form in terms of slightly redefined variables. 
Let 

Ai = a { + ViX 



8 Note also that we have used conventions in which 

[V p , Va]A„ = A RP pa , RP e0 = R aB 
With these conventions the scalar curvature of a unit d sphere is d(d — 1) and the Ricci tensor of a unit sphere is given 
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where \ isi chosen to ensure that V^a* = 0. This is the usual split of a gauge field into its pure curl 
and pure divergence parts. Note that 

fij = Q% &j 9j 0>i Fij . 

We also define the effective pressure 



in terms of which (2.8) reduces to 

V 4 p e + d QVl + v.Vv t - v (V 2 ^ + R ijV j ) = -d a l - v\f ri (2.9) 

(2.9) is precisely the Navier-Stokes equation with forcing function generated by an effective background 
electromagnetic field (with a — and spatial vector a 1 ) on the effectively charged fluid. 

2.3 Charged Fluids 

If the fluid under study carries an extra set of conserved charges, it obeys an extra set of conservation 
equations of the form V M J M = 0. It is easy to verify that these equations all reduce to the condition 
that the fluid velocity is divergence free in the scaling limit under study in this section. Consequently 
charged fluids obey the same equations as uncharged fluids in the scaling limit under study in this 
subsection. 

2.4 Reduction of Cauchy Data 

Returning to the study of uncharged fluids, the Cauchy data 9 of the parent relativistic fluid dynamical 
equations consists of d real functions of space; the value of the pressure field and the value of the d — 1 
independent velocity fields on an initial time slice. As the fluid dynamical equations are of first order 
in time, the Cauchy data of the problem does not include the time derivatives of all these fields. 

Now let us examine the Cauchy data of the incompressible Navier-Stokes equations. Note that 
the spatial divergence of (2.8) 

V 2 p e = -V^VjV - vViZy + V, [(-vR) + /]) v j ] (2.10) 

determines the pressure of the fluid as a function of the fluid velocity (only the velocity - not its time 
derivatives ). In other words the independent data in the fluid is simply given by d — 2 real functions 
that parameterize an arbitrary divergence free velocity field. 

It follows that two of the degrees of freedom of the equations of relativistic fluid dynamics are 
lost on taking the scaling limit of the previous subsection. These two degrees of freedom are simply 
the fluctuations of the pressure and the divergence of the velocity. At the linearized level these two 
degrees of freedom combine together in sound mode fluctuations. (Note that the relativistic dispersion 
of sound implies that a sound mode has twice as much data as a shear mode.). Consequently the 
reduction of Cauchy data in the scaling limit of this paper follows simply a nonlinear restatement of 
the observation that sound waves are projected out in our scaling limit. 

9 We thank S. Trivcdi for very useful discussions on the topic of this subsection. 



Pe=p- 7:^00 + X 
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3. Symmetries 



As we have noted above, the Navier-Stokes equations may be obtained by applying the appropriate 
scalings to the equations of relativistic hydrodynamics with an arbitrary equation of state. The 
parent hydrodynamical system may in particular be chosen to be conformally invariant. It is natural 
to wonder whether either the dilatation or special conformal transformations of the parent theory 
descend to a symmetry of the the Navier-Stokes equations. In this section we answer this question in 
the affirmative. 10 

3.1 Dilatations 

Let us first explain how this works for dilatations. Dilatations consist of a diffcomorphism (x'Y = 
T' = T, (u'Y = ^j- {g')nv = X 2 g llv compounded with the Weyl transformation x^ = (x'Y, g^ v = 
^F, = \{u'Y, f = XT. In sum 

= y , 9^ = «" = u", f = XT. 

This action on coordinates, temperatures and fields, is a symmetry of the equations of conformal 
relativistic fluid dynamics. While this operation commutes with our scaling limit, it is not by itself 
a symmetry of the Navier-Stokes equations. This is because the kinematical viscosity v, regarded as 
a parameter in the Navier-Stokes equations, is proportional to ^ and so changes under the variable 
transformation described above. The Navier-Stokes equations are however simple enough to allow 
dilatation transformations listed above to be modified into a true symmetry of the Navier-Stokes 
equations by 'absorbing' the transformation of v into 'anomalous' transformations of time and velocity. 
The result of this procedure is simply (1.3) with which we started this note. 

3.2 Special Conformal Transformations 

The situation is more straightforward for special conformal transformations. The scaling law for the 
velocity and temperature fields, under a special conformal transformation, may also be obtained by 
compounding a diffcomorphism with the appropriate Weyl transformation. Restricting to infinitesimal 
conformal transformations, we find 11 12 

5x» = -2c.xx» + x 2 <? 

Su" = -2 [x^c u - x u cf] u v - 5x u d u u» (3.2) 
ST = 2c.xT - 5x u d u T 

Now note that special conformal transformations induce an additive shift on the temperature 
fluctuation, ST, proportional to x.cT where T is the temperature of the background. In order that 



10 This section was worked out in collaboration with R. Loganayagam. 

11 In order to verify the covariance of local equations under these symmetry transformations, it is sufficient to omit 
the terms proportional to Sx^d^ but instead to transform all derivatives according to the rule 

8 (dp) = 2 [cpx.d — xpc.d + x.cdp\ 

12 Under this transformation, the shift in a conformal stress tensor is given by 

ST' 11 ' = 2d{c.x)T^ v + 2{x x c» - x"c x )T^ + 2{x x c v - x v c x )T£ - 8x x d x T'" / (3.1) 

Upon accounting also for the shift in the derivatives, it is easy to convince oneselves that the equation of energy 
momentum conservation, for any identically traceless stress tensor, is invariant under special conformal transformations. 
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this shift respect the e 2 scaling of ST we are required to scale c oc e 4 and c» cx e 3 . Imposing this 
scaling and retaining terms only to leading order in the velocity expansion (3.3) reduces to 13 



St = 
Sx l = -t 2 d 

(3.3) 

Sv l = -2c l t + t 2 c,diV l 



ST = 2(-c°t + cV)T + t 2 Cj d 3 T 



It is clear that the symmetry generated by c° acts trivially (it does not act on coordinates or 
velocities, but merely generates a shift, linear in time, of the pressure; more about this below). However 
the symmetries generated by c 1 act nontrivially, and are directly verified to be symmetries of the 
Navier-Stokes equations. 14 



3.3 The Full Symmetry Group 



We have thus discovered that the Navier-Stokes equations enjoy invariance under a conformal sym- 
metry group. The generators of this group are the dilatation D, special conformal symmetries Ki 
Gallilian boosts Bi, the generator of time translations (energy) H, momenta Pi and spatial rotations 
Mij . The action of these generators on velocity fields is given by 



Dv j = (-2td t - x m d m - 1) v 3 

K iV j = - 2tSij + t 2 d l v : > 

BiV 3 =Sa — tdivi 

. 3 (3.4) 
Hv J = — dtv 3 

PiV J — — diV 3 

M ik v 3 =SijV k - Skjv* - (x k d l - x l d k )v 3 



13 As above, in verifying the covariance of equations under these transformation, it is sufficient to omit the terms 



S(d t ) = 2tc i d i 8(di) = 0. 



14 Undcr the transformations listed in (3.3) we have 8p e = = 2c.x so that 8dip e = 2c'. Further S^+v.Vv') = —2c' 
and the viscous term is unchanged. Adding all terms together we have a symmetry of the equations. 
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The commutation relations between these various generators is given by 



[D, Ki] 


= - 3Ki 




[D, Bi] 


= -Bi 




[D,H] 


=2H 




[D, Pi] 


=Pi 




[D, M i:j } 


=0 




[M tJ1 P k ] 


= - 5 ik Pj 


+ SjkPi 


[M l0l K k ] 


= - S ik Kj 


+ Sj k Ki 


[M l3 ,B k ] 


= - SikBj 


+ Sj k Bi 


[Mij,H] 


=0 




[Ki,Pj] 


=0 




[K^Bj] 


=0 




[Ki, H] 


= -2B t 




[H,B 3 ] 


= -Pi 





The symmetry algebra listed in (3.5) may presumably be obtained from the appropriate contraction 
of the parent symmetry algebra SO(d, 2). A related study is currently under progress in [21]. 

In addition to the symmetries listed above, the Navicr-Stokes equations have an infinite dimen- 
sional group of trivial symmetries, under which the pressure is simply shifted by an arbitrary function 
of time. These are symmetries of the equation because only gradients of the pressure enter the Navier- 
Stokes equations, and they are trivial because the pressure is not really an independent variable of the 
Navier-Stokes equations, which may in fact be eliminated by taking the curl of those equations. For 
this reason one need not keep track of the action of symmetry generators on the pressure. However it 
is not difficult to do so and we find 



D Pe =(-2td t ~x m d m -2) Pe 
K iPe =2x l + t 2 d lPe 

BiPe = — tdiV e 

(3.6) 

Hp e = - d t Pe 
PiPe = - dip e 

Mi k p e = - {x h di - x l d k )p e 

This action of the symmetry generators on the pressure field do not quite yield the commutation 
relations (3.5), but instead have additional terms on the RHS corresponding to generators of the 
trivial symmetries referred to above (i.e. generators that shift the pressure field by a function of time). 
Spatial derivatives of the pressure field, however, honestly transform according to the algebra (3.5). 
Consequently, the symmetry algebra (3.5) is not represented on the pressure field itself, but on its 
spatial derivatives. This is vaguely reminiscent of the fact that the two dimensional conformal group 
has well defined action on all derivatives of massless scalar fields, but not on the field itself. 
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4. The gravity solution dual to fluid dynamics in the scaling regime 



4.1 The dual bulk metric 

As we have described above, any solution of the incompressible non-relativistic Navier-Stokes equations 
solves the equations of fluid dynamics dual to gravity up to 0(e 3 ), under the scaling listed in the 
previous section. Now in [9, 10, 13] the equations of fluid dynamics are obtained as the Einstein 
constraint equation of a bulk asymptotically locally AdS^+i space. It is thus natural to define the 
gravitational dual to a solution of the Navier-Stokes equations as any small fluctuation about a black 
brane background that solves all of Einstein's equations (constraint as well as dynamical) to cubic 
order in e. Adopting this definition, it is easy to read off the bulk metric dual to the Navier-Stokes 
equations from an appropriate scaling of the bulk metric of [9, 10, 13]. 

In computing the bulk metric upto C(e 3 ) in the sense described above, it turns out that terms 
from only the zeroth order and the first order in the derivative expansion of the gravitational solutions 
of [9, 10, 13] are relevant. The metric up to first order in derivative has the following form 

ds 2 = dsQ + ds\ 
where 

cLsq = —2u fl dx ti dr + ^ d2 u ll u v dx* 1 dx u + r 2 g fJ/l/ dx tt dx v 



2 

d~ 1' 



ds\ = -2ru v (u a V a ) u^dx^dx" + - — -r (V a u a ) u^ l u v dx lx dx v + 2br 2 F {br) a^dx^dx" 



F(x) = / dy 



d-l _ i 



y(y d -l) 



b=^ f = b» + Sb 
T = T + ST 

Here V denotes the covariant derivative with respect to the full boundary metric which is equal to 
a background g^ v plus perturbation H^ v . T is the temperature of the background blackbrane. The 
terms that appear in the metric involve covariant derivatives of the d velocity with respect to the 
full boundary metric. These terms can be expressed as covariant derivatives of the d — 1 velocity Vi 
and the metric perturbation Ai = H oi with respect to spatial part of the background metric g^ . 

V iUj = ViVj + 0(e 4 ) 

1 1 

ViU + V u 4 = d (vi + Ai) - -dihoo - ^di{vjV 3 ) - v^F^ + C(e 4 ) 

v>" = Vjv j + 0(e 4 ) 
^ v> - 0(e 4 ) 

^ v> 4 = d { Vl + A t ) - ^dihoo + {v j Vj)vi - v j Fij + G(e 4 ) 

Fij = diAj djAi 

Here V denotes the covariant derivative with respect to g^ . The raising and lowering of the i, j 
indices are also with respect to the metric gij. To simplify the expression of cr M „ in (4.2) the constraint 



(4.2) 
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WiV 1 — has been used. Using these expressions the derivative part of the metric can be written as 



ds 2 = b Q r 2 F(b r) (V tVj + VjVi) dx i dx j - 2b r 2 F (b r)v j (V l v + V v t ) dt dx l 
+ 2r (do(vi + A t ) - ^dih 00 - v j Fij + (v j Vj)v^ dt dx l + G(e 4 ) 

Here the first term is of order e 2 and the last two terms are of order e 3 . Since from the constraint 
equations ViV 1 = 0, there is no contribution from the scalar sector. The zeroth order metric can also 
be expanded in powers of e. It turns out that to solve Einstein equation up to order e 3 ,it is sufficient 
to expand the zeroth order metric up to order e 2 in the fluctuations. 

ds l = j \ o dt 2 + f 2 {- dt2 + gindx^x 3 ) + 2dt dr 
2 

(Ai + Vi) dt dx 1 - 2 (Ai + Vi) dx l dr 



ljd r d-2 



+ d+ l d _ 2 (-d 5b + vjv j - /i o) dt 2 + (Ai + Vi) (Aj + Vj) dx l dx j - (-VjV j + h o) dt dr 
b r r 

(4.4) 

Here the first line is of order e°, the second line is of order e 1 and the third is of order e 2 . The full 
metric is given by the sum of (4.4) and (4.3); and solves Einsteins equations to 0(e 3 ) provided the 
velocity and temperature fields above obey the incompressible Navier-Stokes equations (1.1). These 
equations imply in particular that 

^ = -Vii^'vy - uV'iZy + V, [(-vR) + Fj) v 3 ] + ^SJ 2 h m - d (V.A) (4.5) 

an equation that determines ST (and hence Sb in (4.4)) as a spatially nonlocal but temporally ultralocal 
functional of the velocity fields v l . It follows that though the bulk metric at x^ is determined locally as 
a function of temperatures and velocities at x^ , it is not determined locally as a function of velocities 
at x^ . This non locality is a consequence of the infinite speed of sound (and consequent action at a 
distance) in our scaling limit. 

4.2 Comments on the bulk metric 

Several comments on the bulk metric (4.3) and (4.4) are in order. Let us first spell out some termi- 
nology. We refer to the parts of the bulk metric that are proportional to a linear combination of dt 2 
or J^iidx 1 ) 2 as its scalar components. Terms in the metric proportional to dtdx 1 are called its vector 
components, while terms proportional to dx % dx 3 are referred to as its tensor components. 

Now the derivative expansion of [4-13] solves the Einstein equations in a multi step process. In the 
first step the dynamical Einstein equations are used to determine the r dependence of the bulk metric 
as a function of boundary data, in each of the sectors described above. The requirement of regularity 
of the future horizon is then used to determine the boundary data of the tensor sector in terms of the 
boundary data in the vector and scalar sectors (roughly the fluid temperature and velocity) . 

In the scaling limit described in the previous subsection, the tensor sector of the bulk metric is 
particularly simple. It is given by 

b r 2 F(b r) (\7 lVj + \7 jVi ) dx l dx j + — L (A t + Vi ) (Aj + v^dx'dx 3 . (4.6) 
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Note that even though the velocity Vi is scaled to zero in the scaling limit under study in this paper, 
the leading order in e metric (4.6) includes a quadratically nonlinear in the velocities. This is a 
consequence of the fact that we scale distances to infinity at the same rate as velocities to zero. Had 
we simply scaled amplitudes to be small, while keeping distance scales finite, we would have ended up 
with only the first of these two terms in (4.6) and the resultant geometry would simply have been the 
dual of linearized fluid dynamics. 

The Navier-Stokes equations are obtained out of the Einstein constraint equations acting on the 
bulk metric described above. The first term in (4.6) gives rise to the viscous term in the Navier- 
Stokes equations, while the second term in (4.6) is the origin of the nonlinear convective term in 
those equations. It is consequently not surprising that ratio of the first to the second term in (4.6) is 
proportional to the Reynolds number Re of the flow: in fact it is of order Re x r d F i r \ ■ Now F(r) ~ -. 
at large r; so that the viscous term in (4.6) appears to dominates the nonlinear term in that expression 
when r > (i?e) 1//d_1 . This appearance is atleast partly a coordinate artifact, as is evidenced by the 
fact that the two terms in (4.6) contribute equally to the Einstein constraint equations at every r (this 
follows as as these equations are independent of r). 15 

5. Simple Steady State Shear flows at arbitrary Reynolds Number 

Consider a fluid on J?** -1 ' 1 , subject to an effective forcing 

a x = dujdka(k,uS) exp[iky + iu)t] 

J (5.1) 

a t = 0, (i^ x) 

Here x and y parameterize orthogonal displacements in the d—1 spatial directions, and a(k, ui) is any 
function. 

In the presence of this forcing function, it is easy to verify that the velocity field 



v x = - dkdui- 

J (5.2) 

Vi = i ^ x 

(together with the pressure field which may be obtained by integrating (2.10) ), is an exact solution 
of (2.9). 

In order to be concrete let us make the simple choice a(k, uS) = a(—k, —u>) = aS(k — ko)S(ui — luq). 
In this special case the forcing function is given by 

a x — a cxp [ik y + iuo^t] + cc 
Oj =0, (i^ x) 



(5.3) 



and the velocity field is given by 



ae ik y+iuot 

1-^ (5-4) 



Vi = i ^ x 



15 The ratio computed above naively suggests that nonlincaritics at different effective Reynolds numbers are separated 
in the bulk radial coordinate; a result reminiscent of the approximate locality in scale of high Reynolds number fluid 
flows. While this interpretation is too quick for the purposes of evaluating the Einstein constraint equation, perhaps 
the naive ratio of the two terms in (4.6) is physically relevant for some physically interesting questions. We feel this 
deserves further contemplation. 
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16 This solution is characterized by two dimensionless numbers, a — j^l. an d 5 = £rL_ The Reynolds 
number for this solution is given by 

and can be made large by taking a to infinity at fixed b. This may be achieved, for instance, by 
increasing the strength of the forcing amplitude a keeping everything else fixed. 

In Appendix A we outline a linear stability analysis of the special flow described above. We find 
that in d = 2 this flow is unstable at large enough Reynolds numbers for a range of the parameter 
b. As we have described in the introduction, this suggests that the flow is dynamically interesting, 
and likely turbulent, at asymptotically high Reynolds numbers. Assuming this is the case, we have 
identified atleast one gravitational system that is dual to steady state turbulence. In some detail, 
consider gravity with a negative cosmological constant, subject to the small z boundary condition on 
the metric 

lim ds 2 = (dz 2 — dt 2 + edtdx l a,i + dx 1 dxi) 

z^O Z l y ' 

(where i = i . . . d — 1 and a, is the gauge field (5.1)). The steady state solution of this gravitational 
system is dual to a turbulent flow of d dimensional fluid dynamics whenever the parameters in the 
function dj above are chosen to satisfy the inequality a » \f\ + b 2 . 

In Appendix B we have presented a generalization of the solution described in this section to the 
forced flow of a fluid on a sphere. In that situation one should obtain the dual of a turbulent fluid in 
an spacetime that is asymptotically global AdS. 
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Appendices 

A. Stability Analysis 

In this appendix we study linear fluctuations about the solution (5.4) for the special case d = 3. In 
two spatial dimensions the dual of a divergenceless velocity is the gradient of a scalar, and the field 
strength of a gauge field is itself dual to a scalar. In other words we can write 

ViyAx (A1) 

fij — Cijf 



16 This special solution may be thought of as living either on a spatial R d 1 or on a spatial torus, whose y edge is of 
2tt 



size t— 
fe 



- 15 - 



The Navier-Stokes equations (2.9) may be rewritten (after eliminating the pressure by taking the 
curl) in terms of these scalars as 

-V 2 x + ^V 4 X - e- ?m V mX V J V 2 X + / + e^V tX V q f = (A.2) 

In terms of these variables the solution (5.4) takes the form 

fo = — ifcoaexp [ikoy + itoot] + cc 

1 ae ik y+ivot r A3 \ 

Xo = -— x + cc 

1 — 



A.l Setting up the Eigenvalue Equations 

Let us now set 



X = Xo + ' 



(e iqax j du>dk Xqo {cj,k)e iu ' y+ikx +e- iqox J divdk Xqo {cu,k)e- iojy - ikx 



where 9 is a small parameter, and solve (A.2) to linear order in 6. Note that while our fluctuation has 
a specific frequency (namely go) m the x direction, it is a sum over frequencies in the y direction and 
in time. This is necessary as our background breaks translational invariance in the y and t directions, 
but preserves this invariance in the x direction. It is easy to work out the linear equations that our 
fluctuation coefficients Xqo {u>,k) obey. We find 

[iu(k 2 + ql) + v(k 2 + qlf] X<a,(w,«) + [/3{iq ) ((fc - k ) 2 + q 2 - k 2 ) + aiq k 2 ] Xq „(u -Uo,k- k ) 
+ [f3* {iq ) ({k + k ) 2 + q 2 - k 2 ) + a*iq k 2 ] Xqo (w + u , k + k ) = 

(A.4) 

where 



Let us define 

fk(l, k) = Xqa [w (7 + k )i fc o(« + k)] . 

The fluctuation equation (A.4) may be recast in terms of fk as (we will usually omit to write the 
functional dependence of fk in the equations that follow) 



+ k)((K + k) 2 + c 2 ) + ((k + k) 2 + c 2 ) 2 

. , ( (K + k+l) 2 +C 2 -l , \ 



, , • | (K + k-l) 2 + C 2 -l 

fk+ica [ — — h 1 ) fk-i 



(A.5) 



The dimensionless quantities a, b and c are parameters in (A.5) (a and b were defined in section 4 
while c = |^). 7 and k are variables in this equation; (A.5) is the condition that a matrix acting on 
the columns {fk} has a zero eigenvalue. Together with the physical requirement that fk decay at large 
|fc|; this eigenvalue equation yields an expression for unknown temporal frequency 7 as a function of 
k. 17 That is, the equation (A.5) should yield a dispersion relation of the form 

7 = i(n,a,b,c,ri) (A.6) 



17 Recall that k is a real number (in order that the mode in question is well defined at all y) of unit periodicity. 
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where the integer n labels which of the infinitely many solutions to the zero eigenvalue condition we 
have chosen to study. A result for 7 with a negative imaginary part represents an instability of the 
system. 

Let us study the large |fc| asymptotics of the variables fk in a little more detail. According to 
our boundary conditions, at large positive k fk+i is much smaller than fk-i- As a consequence (A. 5) 
implies that 

h - ica * *r 2( r^ *fi+o(A)) (a.7) 



fk-i l-ib k 2 (k + 4n + i) \ y k 2 
It follows that, at large k, 

n , ( ica \ k Y(k + 2k - 1) 

fk ~ D(a,b,CK) r — ^ -. (A.8) 

J v '' ; V*- i& / r(fc + i) 2 r(fc + 4 K + f + 1) v ; 

where D is a constant. This estimate is valid provided that k 3> 1 and that /c 2 ^> ^=p ■ Similarly, 
the behavior of /fe at large negative values of fc is given by 

f fe » D'(a* , 6, c, «) f " r( - fe - 2l6 - 1} ■ (A.9) 

J v y V* + *by r(-/c + i) 2 r(-fc-4 K + i + i) v ; 

where D' is an independent constant. Note \ is re& l provided that D(a, b, c, n) = D'(a* ,b, —c, —k); a 
condition that we consequently demand on physical grounds. 

We are interested in determining the dispersion relation (A. 6) as the condition for the existence of 
a solution of (A. 5) that achieves both asymptotic conditions (A.8) and (A.9) (this leads to an equation 
as, generic solutions of (A. 5) that obey (A.8) blow up at large negative k). 

A. 2 Qualitative nature of the solution 

Let us first note that the off diagonal terms in (A. 5) are both proportional to ax c. The proportionality 
to a reflects the fact that the problem gets strongly coupled at high Reynolds numbers; we will have 
a lot more to say about this below. The fact that these terms are proportional to c, however, implies 
that fluctuations with no x momentum are governed by the linear Navier-Stokes equations, and so 
are always stable. An instability, if it occurs, must necessarily break a symmetry (in this case of x 
translations) that the background solution preserves. 

In much of the rest of this section we specialize our analysis to n = and c = 1 (i.e. ko = qo) but 
for a reasonably wide range of the parameters a and b. We perform this specialization for convenience 
in the numerical analysis that we will describe below; we expect qualitatively similar results at all 
values of k and reasonable nonzero values of c, though we have not explicitly tested this expectation. 

Let us start by painting a qualitative picture of physically relevant solutions of (A. 5). We start with 
the simple limit a — > 0. Solutions are given exactly by fk = 5k.K with 7 = — K + ib(K 2 + 1) for K rang- 
ing over integers. The time dependence of this solution is given by /k(t) = fk(0)e~ lKuJot ~ buJo ( K + 1 )*; 
it follows that all these modes decay in time, so that our solution is stable against small fluctuations. 

At small but nonzero values of Reynolds number Re = ^-^ +b2 the solutions to (A. 5) are small 
perturbations of the solutions described in the paragraph above. On the K th such solution fk is 
nonzero at k = K, but decays rapidly as k moves away from K. In particular fk-K oc (Re)\ k ~ K \ . 18 
The dispersion relation for 7, described in the previous paragraph, is corrected in a power series in 
(Re) 2 . We will compute the first term in this series below. 



s This decay is visible, for instance, in the asymptotic forms (A.8) and (A.9). 
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As the Reynolds number is increased, the width in the distribution of ff. (as a function of k 
centered around k = K) increases. For Re 3> 1 and Re 3> \pn this width is of order k ~ (i?e)5. (see 
(A.8) and (A.9)). 

In summary, the fluctuation mode we study is highly localized about a particular y momentum 
at small Reynolds number, but consists of a 'cascade' or roughly equal superpositions of order (Re)^ 
modes at large Reynolds numbers. The fluctuation modes oscillate rather than growing exponentially 
in time at small Reynolds numbers; however the stability properties of these modes must be elucidated 
by explicit calculation at large Reynolds numbers. 

A. 3 Numerical Evaluation of the frequency 

In this subsection we present the results of our rough numerical evaluation of the frequency 7 of the 
mode with n = 1 as a function of the Reynolds number, at various values of b. The method we use 
is very simple. We start with a value of k large enough for the asymptotic form (A.8) to be valid 19 
. We then use the recursion relations (A. 5) to evaluate /i//o- We then independently evaluate the 
same ratio using (A.9) and recursion relations. Equating these expressions for the same ratio gives us 
an equation which we use to solve for 7 as a function of all other parameters. Of course this equation 
has several solutions. At small Reynolds number we choose the solution that is near to 7 = — 1 + 2ib 
(i.e. the mode with K = 1 in the language of the previous subsection) , and then 'follow' this root as 
we numerically increase the Reynolds number (in small steps) to large values. We have performed all 
our calculations using Mathematica. 



Imaginary part of y 
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Figure 1: The imaginary part of the frequency 7 plotted against a, for k = 0, c = 1 and the root with K = 1. 
Note that 7771(7) = 2 at b — and decreases monotonically as |6| is increased. 

Let us describe our results in detail at b = 1. In Fig. 1. below we present a plot of the 
imaginary part of 7 as a function of Reynolds number. In Fig. 2 below we plot the real part of the 
same frequency as a a function of Reynolds number. Notice that the imaginary part changes sign 
(indicating a transition to instability) at a Reynolds number of order 13. 

We have performed the same analysis at several different values of b. At each value of b the 
imaginary part of 7 turns negative at a particular value of the Reynolds number. In Fig. 3 below we 
have plotted this critical Reynolds number as a function of b for a range of values of b. 

19 We have performed computations with values of this starting k ranging from 20 to 8. The plots presented below are 
generated using the starting value 10. We have verified that our starting value of k is large enough, by checking that 
our results are not significantly affected by increasing k. 
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Real part of y 
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Figure 2: The real part of the frequency 7 plotted against a, for k = 0, c = 1 and the root with K = 1. Note 
that Re(-/) = —1 at 6 = and decreases monotonically as |6| is increased. 



Re = a/Sqil[l+b A 2] 




Figure 3: The curve that separates the region of stability (below) from instability (above) of the eigenvalue 
at K = 1, plotted as a function of the Reynolds number on the y axis versus b on the x axis. We have used 
k = 0, c = 1 in the calculations that generated this plot. 



As our principle aim in this subsection is to establish that the small fluctuations about our solutions 
are unstable at high enough Reynolds numbers we have not attempted to carefully estimate the errors 
in our numerical calculations. However we think it is unlikely that the errors in, for instance, Fig. 3 
exceed a few percent. 

As we have described above, at every value of b we start our calculations at small Reynolds 
numbers which are then slowly increased. As a check on our numerics, at small we have compared our 
numerical results for 7 versus k against the predictions of perturbation theory (see the next subsection 
for details) at lowest order in (Re) 2 ; we find good agreement. In fact at b = 1 and a — 1/100, 
the difference between our numerical result and the zeroth order answer matched the prediction of 
perturbation theory upto one part in 10 5 ; a more accurate result than we had expected. 
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A. 4 Perturbation theory at small Reynolds numbers 

Let us formally identify the equation (A. 5) with free index k with 

(k\M\tp) = ^2(k\M\m)(m\ip) 



and simultaneously identify 

(«#) - fk 

With these formal identifications, the set of equations (A. 5) (for all k) are simply the equation for the 
operator M to have a zero eigenvalue. The corresponding eigenvector is specified by the values of f k 
on this solution. 

Let us write M = Mo + Mi where Mo is a diagonal operator in the basis described above and Mi 
is an off-diagonal operator, proportional to a. The elements of the matrix M and Mi are given by 



(fk\M \f k )=A(k) 
i 



2\2 



( 7 + fc)((« + k) 1 + c 2 ) + ((« + k) 2 + c 2 ) 



{f k \M 1 \f k _ l )=B{k) 

. ( ( K + fc-l) 2 +c 2 -l 

{f k \M 1 \f k+1 ) = C{k)5 Uj+l) 

( K + k + l) 2 +c 2 - 1 



(A.10) 



l + ib 



+ 1*. 



3fe,(j+l) 



When a = the eigenvectors of the operator M are simply \K) with eigenvalue A(i^). This 
eigenvalue vanishes when 

7 o - -K + i b [(k + X) 2 + c 2 ] . 

We now wish to study how this eigenvector - and the corresponding solution for 7 that keeps the 
eigenvalue zero - evolves in perturbation theory in a. 

To lowest nontrivial order in perturbation theory, the shift in the eigenvalue A(k) is given by 



A{k) = A(k,u) + (fk\ M i\h-i)(fk-i\Mi\f k ) + (fk\Mi\f k+ i) (fk+i\Mi\f k ) 



A(k)-A(k-1) 



A(k)-A(k + 1) 



(A.11) 



A(k)-A(k-1) A(k)-A(k + 1) 
In order that the eigenvalue vanish at this order we must have 

i b 



1 



-K + i b [(k + K ) 2 + c 2 ] 



(K + k) 2 + c 2 



B(K)C{K - 1) C{K)B{K + 1) 



A(K-l) 



A(K + 1) 



where the third term is evaluated at 7 = 70 and we have used the fact that A{K) = at 7 = 70. 



B. A simple flow on a sphere 

It is not difficult to generalize the simple laminar flow presented in section 4 above to a shear flow 
on a 2 sphere. The corresponding dual gravitational solution to this flow is asymptotic to (slightly 
perturbed) global AdS space. 
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For concreteness we choose d = 3 and choose the base metric of our space to be the two sphere, 

i.e. 

ds 2 = -dt 2 + dO 2 + sin 2 6dct> 2 (B . 1 ) 

Our velocity and forcing function fields will be vector fields on the two sphere, so be briefly pause 
to recall some necessary definitions. Recall that the spherical harmonics, form a basis for the 
expansion of an arbitrary scalar field on the sphere. On the other hand an arbitrary vector field on the 
sphere is given by linear combinations of diY^ and ejdjY^ where the e symbol includes relevant factors 
of ^/g. Now, in our problem, both v % and a 1 are divergenceless. It follows that each of these fields 
may be expanded in a sum over only the vector (rather than also the derivative of scalar) spherical 
harmonics. 

In order to obtain one solution (that can easily be generalized in many ways) to the Navier-Stokes 
equations, let 

ai = ae?d m Yi{e,4>). 

The velocity configuration 20 

"* = - - .,° + i)- a xeWftM). 
1 — i K ' — 

(together with the implied pressure field) yields a steady state, laminar, shear flow. The Reynolds 
number of this flow is given by (5.5) with a and b now given by a = and b = t/ ('('+ 1 )~ 2 _ As for the 
solution presented in section 4, we expect this flow to be unstable to small perturbations at fixed b in 
the limit of large a (we have not performed the linear stability analysis in this case). 
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